Quantification of ambient PM2.5 concentrations adjacent to informal brick kilns in the Vhembe District using low-cost sensors

The widespread exposure to ambient PM2.5 poses a substantial health risk globally, with a more pronounced impact on low- to medium-income nations. This study investigates the spatiotemporal distribution of PM2.5 in the communities hosting informal brickmaking industries in Vhembe District. Utilizing Dylos DC1700, continuous monitoring of PM2.5 was conducted at nine stations adjacent to informal brick kilns from March 2021 to February 2022. The study determined the correction factor for PM2.5 measurements obtained from the Dylos DC1700 when it was collocated with the GRIMM Environmental Dust Monitor 180. Additionally, the diurnal and seasonal variations across monitoring stations were assessed, and potential PM2.5 sources were identified. The study also evaluated the compliance of ambient PM2.5 concentrations across the stations with the South African National Ambient Air Quality Standard (NAAQS) limits. Annual PM2.5 concentrations for the stations ranged from 22.6 to 36.2 μgm−3. Diurnal patterns exhibited peak concentrations in the morning and evening, while seasonal variations showed higher concentrations in winter and lower concentrations in summer and spring. All monitoring stations reported the highest daily exceedance with respect to the daily NAAQS limit in the winter. Major PM2.5 sources included domestic biomass combustion, vehicular emissions, industrial emissions, and construction sites. Well-calibrated low-cost sensors could be employed in suburb regions with scarce air quality data. Findings from the study could be used for developing mitigation strategies to reduce health risks associated with PM2.5 exposure in the area.

www.nature.com/scientificreports/Attempts to collect ambient real-time air quality data with high-resolution is often encumbered by the limitations of the air quality instruments 21 .Ambient air quality monitoring is often assessed using a limited number of stationary air quality monitoring stations, basically because the widespread deployment of these regulatorygrade air monitoring sensors is too expensive 1, [22][23][24] .In the developing world, especially in Africa, the use of the conventional methods of air quality monitoring is too sparse to meet the growing challenges of development 25,26 .This toughens assessing the gravity of the air pollution problem, as there is limited or no evidence to show that air quality is a challenge, thus creating a harmful cycle 23,27 .In South Africa, the available air monitoring stations are mostly located in and around major industrial and densely populated urban centers 25,28,29 .Several studies in South Africa have also assessed ambient PM 2.5 concentrations using a single or network of gravimetric equipment around low-income settlements [30][31][32] .The Vhembe District of Limpopo, on the other hand, do not receive much preference.Although, particulate matter monitoring occurs in some districts in Limpopo, the Vhembe District currently has no functional ambient air quality monitoring station for PM 2.5 29 .To generate a more robust PM 2.5 database with high temporal-spatial resolution, a relatively inexpensive and readily operated real-time measuring apparatus for PM 2.5 monitoring is essential.
By addressing cost concerns and incorporating portability, the low-cost sensors (LCS) are increasingly gaining global acceptance in their use for air quality assessment [33][34][35][36] .The easy operation of LCS promotes public participation 37 , easy accessibility to air quality data and public awareness 38 .Additionally, LCS enable extensive monitoring of air pollution with higher spatiotemporal resolution, thus augmenting data gathered through conventional technology 22,39 .The robust PM 2.5 database generated using PM LCS, enables policymakers to target interventions at highly polluted regions, ultimately improving air quality 40 .Dylos DC1700 (Dylos Corporation, Riverside, CA) is one of the commonly used LCS for monitoring particulate matter.Dylos being an optical particle counter operates by illuminating a single particle at a time to count the number of particles in defined size bins 41 .The primary challenge with light-scattering instruments is that the instrument response relies upon both the size distribution and the number of particles, instead of the cumulative mass of aerosol.This limitation could be minimised, and reliability enhanced through periodic calibration against a federal reference method 36,42,43 .Dylos is also limited by its detection limit of 5 µm, which makes it impossible to measure ultrafine particles that pose a greater risk to humans 34,44 .Dylos is most conveniently used for indoor air quality monitoring due to the short battery life and its fragility to withstand outdoor environment.The use of Dylos for monitoring indoor air quality has been reported in a number of studies 22,34,45,46 .When used for outdoor air quality monitoring, the study is often conducted in partially enclosed environments, or used in the outdoor environment for relatively short period of time during dry weather [47][48][49] .To guarantee the sensor's long-term effectiveness, further investigation on Dylos stability is required 38 .
Brickmaking in many developing nations is acknowledged as a polluting industry characterised by small-scale enterprises that are challenging to monitor and regulate [50][51][52][53][54] .In 2010, the clamp kiln, a commonly used brickmaking technology was enlisted among the activities requiring an Atmospheric Emissions Licensing (AEL) for operation in South Africa 55 .The brickmaking industry is known to be one of the thriving commercial activities in the Vhembe District.However, there is limited knowledge on the impact this industry has on ambient PM 2.5 in communities close to the brick kilns.This study is born out of the limited database on PM 2.5 level in ambient air around the brickmaking industry in Vhembe District.The study assessed the ambient PM 2.5 around the informal clamp kilns, using a network of Dylos DC1700 PM sensor for a period of 1 year.For the first time in Africa, a correction factor for PM 2.5 readings from Dylos DC1700 PM sensor was generated.The research also examined the extended suitability of Dylos DC1700 for outdoor monitoring of ambient PM 2.5 concentrations.
The objectives of the study include the following: • To conduct 1-year continuous monitoring of PM 2.5 around the informal brick kilns in the Vhembe District from March 2021 to February 2022, using a network of well calibrated Dylos DC1700 PM monitors.• To examine the compliance of the PM 2.5 levels in the host communities of the brick kilns to the South African National Ambient Air Quality Standard (NAAQS).• To determine the temporal variation of PM 2.5 level.
• To investigate other potential sources of the PM 2.5 in the study area.

Methodology The study area
The study area for this work is the Vhembe District in northern Limpopo province (22.7696°S, 29.9741°E),South Africa as shown in Fig. 1.Through the Kruger National Park, Vhembe District shares borders with Zimbabwe and Botswana in the north-west and Mozambique in the south-east.The area is separated from its bordering countries by the Limpopo River.The district is divided into four municipalities: Mahkado, Collins Chabane, Musina, and Thulamela.It has a vast amount of territory that spans 25,596 km 2 .Its population is estimated at 1.5 million people.The populace primarily relies on mining and agriculture for its subsistence 56 .Vhembe District receives most of its rainfall between October and March 57 .The district receives between 400 and 1100 mm of precipitation annually.The district also experiences the highest temperature during the wet season 57 .Vhembe District has many informal bricks making industries, employing the use of the clamp kiln for firing of bricks.The brickmaking industries are often in clusters, and they are mainly situated in Mannini and Tshilungoma villages.These areas of the district are desirable places for brick producers due to the presence of clay deposits and water bodies 51,58 .

Site selection and instrument installation
Before the commencement of the study, Ethical clearance (SES/19/ERM/08/1309) was obtained from the University of Venda Research and Ethics Committee.Nine locations within 2 km radius from the clusters of informal brick kilns were selected for installation of PM 2.5 monitoring equipment.The monitoring sites were selected based on approval from the municipality, ward counsellors, traditional rulers, and other constituted authorities/ landowners.Another important factor considered for selecting the monitoring sites was the security of the monitoring instruments to be installed.The location of the PM 2.5 monitoring stations were selected to span a range of land use characteristics, including governmental facilities/organizations, academic environments, recreation centre, site near local and arterial roads, and residential areas (Table 1).
Prior to Dylos DC1700 (Dylos Corporation, Riverside, CA) installation, the Dylos units were modified by the technical crew of the Climatology Research Group at the North-West University, Potchefstroom, South Africa.A solar panel, solar controller, higher voltage battery (12V-12A), and a waterproof housing was added to each of the Dylos devices to enhance their performance in the outdoor environment (Fig. 2a).One of the modified Dylos units was then collocated with a GRIMM Environmental Dust Monitor 180 at the North-West University for a period of 5 months.A correction factor was generated from the collocation exercise of the Dylos device  with the GRIMM instrument, and this factor was used to estimate the Particle Mass Concentration (PMC) in µgm −3 in the study.
For the installation of the monitoring devices at the selected points, EPA Probe and Monitoring Path Siting Criteria for Ambient Air Quality Monitoring was followed as practicable as possible.The instruments were installed at 1.5 m above ground level (Fig. 2b), about 1 m from wall/fence and 10 m from the drip line of trees 59 .Height of 1.5 m above ground level is assumed to be the statistical average breathing height of pedestrians [60][61][62][63] .One-year continuous monitoring of PM 2.5 was conducted from March 2021 to February 2022.For the study, hourly meteorological measurements for the study period (March 2021-February 2022) were obtained from the South African Weather Stations (SAWS) situated within the study area (Latitude − 22.9160 and Longitude 30.4050).The meteorological parameters used for the study are the wind direction and wind speed.

Dylos DC1700 instrumentation
The Dylos device is programmed to collect air quality data from the sampled air stream every minute.Data were collected and stored in the memory throughout the stipulated monitoring duration.Due to the short battery capacity of Dylos DC1700, the internal battery of the device was complemented by a solar system with 12 voltage battery.Readings on the Dylos device are displayed in 2 size bins; one counting all particles greater than 0.5 µm (the detection limit) and the other counting only coarse particles above 2.5 µm 34 .Particle number count (PNC) between 0.5 and 2.5 µm is calculated from the difference between the PNC of the large (> 2.5 µm) and small (> 0.5 µm) particles 22,34 .Measurement with this device is expressed in particle count per 0.01 cubic foot (0.283 L) of air 64 .Given that there is basically no epidemiological proof connecting PNC with health impacts, and to communicate the implications of exposure to poor air quality to society 64,65 , it is essential to report the measurement from the Dylos device in the corresponding mass concentration.The corresponding mass concentration could be estimated using the correction factor generated during the collocation experiment 34,66 .The device is provided with exclusive programming Dylos Logger (V.1.60)to empower the logged information from the instrument to be downloaded 64 .The measured data are automatically stored in internal memory and downloaded through a serial port to a computer using a 9-pin serial/Universal Serial Bus (USB) cable 64,67 .Due to the short memory capacity of the device, stored data were downloaded manually at least once a week.

Quality control and data analysis
To ensure data quality, the air quality monitoring exercise was subjected to recognized procedures so that the resultant data are representative and comparable.The method described by Aniyikaiye et al. 66 was used.A quality control check was carried out on the collected data by removing the negative and zero readings.The minute-byminute air quality data from the installed Dylos devices were collated and converted to the PMC using the correction factor generated from the collocation exercise.The correction factor was generated using Eqs.(1) and (2).
where PNC PM2.5 is the particle number between 0.5 and 2.5 µm.PNC small is the particle number count > 0.5 µm particles generated by Dylos DC1700.PNC large is the particle number count > 2.5 µm particles generated by Dylos DC1700.PMC is the calculated Particle mass concentration in µgm −3 .
Generally, a dataset comprising of at least 80% of data capture is required for quality assurance and to qualify for data manipulation and summary 68,69 .However, this was lowered to accommodate the modest functionality of the low-cost instrument.For collation of the daily and monthly data, 70% data completeness was used.In

Local/ regional transport analysis of PM 2.5
Meteorology greatly influences the level of pollutants in ambient air.Among the numerous meteorological factors, wind speed and direction play significant role in controlling the pollutants concentrations in the atmosphere.Bivariate polar plots are often employed in illustrating the variability of pollutant concentration with wind speed and wind direction in polar coordinates 70 .In a bivariate polar plot, wind speed is represented as the distance from the origin and wind direction as the angle from the origin.The varied colours in the polar plot represent the average concentrations of the pollutant variable 71,72 .For the study, the local and regional transport analysis of PM 2.5 was determined using bivariate polar plots (BPP).The PM 2.5 hourly data obtained from the monitoring stations as well as the corresponding hourly meteorological data for the period of study were used for the BPP.

Collocation of GRIMM EDM 180 Against Dylos DC1700
A coefficient of determination, R 2 = 0.75 was attained between the 10 min readings and the corresponding 10-min averaged PNC generated by the GRIMM reference instrument and the Dylos device, respectively, during the 5-months collocation exercise held at the North-West university monitoring station in 2021 (Fig. 3).In addition, high reproducibility levels with R 2 values between 0.99 and 1 were attained between the 1-h PM 2.5 average concentrations generated by all Dylos sensors used during a 5-days calibration exercise held at the North-West University in 2020 73 .High R 2 value of pollutant concentrations among network of similar instruments are indication of consistency and strong similarity in the readings from the instruments 74 .With reference to the Dylos PM sensors, the high R 2 value implies that the Dylos units can generate similar set of readings when subjected to the same environmental condition.For this study, the linear equation generated from the collocation exercise between the Grimm EDM 180 and Dylos DC1700 is referred to as the correction factor.The correction factor was determined using Eqs.(1) and ( 2) above.The correction factor was used for the conversion of the PNC from the PM 2.5 monitoring stations to PMC (particle mass concentration) in microgram per cubic meter.

Data
The summary of data availability for all the PM 2.5 monitoring stations within the period of study are presented in Table 2. 70% data completeness is the threshold, data availability with percentages less than 70 are emboldened, while the empty spaces represent results that were not available due to theft of solar panel.For data analysis, the months with less than 70% data were disregarded and discarded.It is worthy to note that no data was available in October 2021 in all the monitoring stations, which was as a result of high rainfall and low sunlight levels.In addition, the low resolution of data collected between August and November was due to battery failure, which was latter replaced resulting in a much-improved dataset in the subsequent months.Based on findings from the study, Dylos performed best when there was no rain, and the batteries are fully charge.Other reasons for gaps in data for the sampling duration (March 2021-February 2022) include; Dylos turning off automatically, inadequate power supply from the solar panel to the equipment due to insufficient sunlight especially during winter (June, July and August) and spring (September, October and November (SON)), dead external battery, malfunctioning instrument, where the instrument was removed for repair, zero readings from the instrument and theft of the instrument.
Based on the prevailing wind direction, the monitoring exercise was conducted at "upwind" and "downwind" locations to the clusters of informal brick kilns.The upwind/downwind monitoring was intended to provide a measure of the contribution of the brick kiln operations to local PM 2.5 concentrations.Therefore, two and seven monitoring sites situated upwind and downwind the brick kilns, respectively, were used to assess the PM 2.5 levels of the host communities of the informal brick kilns.
For easy analysis, the monitoring stations were subdivided into two groups namely: The group upwind the brick making industries (DR and MG).
The group downwind the brick making industries (BM, LB, MM, MR, MT, ST, TP).
Although data quality control considerably decreased the amount of data available for analysis of most of the locations, locations having at least one valid month in a season was assumed to be representative of that season.For monitoring stations situated downwind the informal brick kilns, PM 2.5 concentrations were recorded at all the stations in autumn (March, April and May (MAM)).Conversely, there was low level (< 70%) of data capture at monitoring sites TP, ST, MT and MR in winter, BM, MT and ST in spring as well as BM and LB in summer (December, January and February (DJF)) due to equipment malfunction and theft.

Daily variability pattern of PM 2.5 in Vhembe District
The daily mean mass concentration of PM 2.5 ranges from 18.6 to 28.4, and 23.5 to 31.7 μgm −3 for the groups of monitoring stations situated upwind and downwind the brick kilns, respectively (Table 3).The rows of data italicised in Table 3 are the seasons which were not represented in overall seasonal plots for the monitoring stations due to limited available monthly data because the season do not have at least 1 month with 70% data completeness.The daily average ambient mass concentrations of PM 2.5 at all the monitoring station reached their maximum during winter (June, July and August (JJA)).The concentrations drop to their minimum during summer for majority of the stations, and in spring at stations ST and LB.Maximal concentrations of PM 2.5 were observed across all the monitoring stations in winter due to the maximal operation of the brick kilns, domestic biomass combustion for cooking and space heating and stable atmospheric condition (Table 4).On the other hand, the reduced levels of PM 2.5 experienced during summer and spring could be due to rainfall.In addition, high variability in PM 2.5 concentration were observed in majority of the stations in winter.The relatively high variation in PM 2.5 in winter compared to other seasons could be partly attributed to the occurrence of more severe PM 2.5 pollution events coupled with the regular pollution events common in other seasons.The daily mean mass concentrations of PM 2.5 were compared with the daily NAAQS of 40 μgm −3 for PM 2.5 .Based on the available data, maximum number of daily exceedances to the daily NAAQS limit for PM 2.5 occurred in all the monitoring stations during winter while PM 2.5 levels fell below the NAAQS for PM 2.5 throughout the study area in summer.

Monthly/ seasonal characteristics of PM 2.5
Considering the monthly variation in PM 2.5 in both the upwind and downwind direction of the brick kilns, relatively higher level of PM 2.5 were observed from late autumn (May) and during winter which are periods of expected high pollution.Studies have shown relatively high levels of ambient pollution during winter compared to pollution level in rainy season [75][76][77] .The abrupt increment of wintertime pollution especially PMs concentrations, is primarily ascribed to extensive use of fuel during winter for large-scale heating 31,32,68,78,79 .Previous findings by Aniyikaiye et al. 51 (2021a) in the study area also indicated highest levels of the brickmaking industrial operations during winter.June and July are the key winter months in the study area.PM 2.5 generating activities in this area such as brick firing, domestic biomass combustion for cooking and space heating are highest during autumn and winter.Maximum PM 2.5 level was observed in June in majority of the stations with winter data except LB which occurred in July (Table 4).
Based on the available data, PM 2.5 levels increased from summer < spring < autumn < winter at all the monitoring stations around the brick kilns.During the wet season (summer and spring), brickmaking activities are reduced in Vhembe District due to unconducive weather condition for brick production, and water lodging of the brickmaking areas 51 .Mazumdar et al. 80 also claimed water lodging of quarried land as a major challenge on brick making industries during the rainy season at Palasbari Revenue Circle, India.The reduced brickmaking activity during the wet season could have contributed to the low level of PM 2.5 recorded in summer and spring.www.nature.com/scientificreports/Additionally, the occurrence of favourable meteorological conditions, such as increased temperature and wind speed as well as high precipitation levels during these seasons aid pollutant dispersion and reduction.According to Yuan et al. 81 , precipitation is one reason for low pollutant contamination in the wet season as the pollutants are washed out of the atmosphere.Figure 4a,b present the seasonal PM 2.5 levels for the respective groups.Comparing this study with similar studies in Africa, low to high levels of PM 2.5 were reported (Table 5).There is disparity in air quality results between study conducted by Novela et al. 82 and the current study despite the location of the study areas in the same district municipality.Unlike this study, which was conducted in a

Spatial distribution of PM 2.5 in Vhembe District
The annual distribution of PM 2.5 across the monitoring sites for the study period is presented in Fig. 4c.This study showed that ambient PM 2.5 concentrations range of 22.6 (MT)-36.2µgm −3 (BM).With respect to the compliance levels of the air quality from the monitoring stations to the annual NAAQS of PM 2.5, none of the annual average PM 2.5 concentrations from all the monitoring sites fell below the annual NAAQS level of 20 µgm −3 .DR and MG monitoring stations situated upwind the brick kilns are expected to have the least levels of PM 2.5 , however, annual average concentration of PM 2.5 recorded in these monitoring stations were relatively higher than the annual averages of some monitoring sites (MR, ST, TP, and MT) situated downwind the brick kilns.The relatively higher PM 2.5 concentrations recorded in MG and DR stations could possibly be attributed to the presence of other prominent sources of PM 2.5 (other than the brickmaking industries) in the upwind direction of MG and DR monitoring stations.
Result from the study also showed that median values of majority of the monitoring stations namely; MG, MR, ST, TP, and MT were lower than or approximately equal to the annual threshold limit of 20 µgm −3 .This implies that about 50% of the PM 2.5 levels reported at the listed stations fell below the annual NAAQS threshold limit for PM 2.5 throughout the monitoring period.On the other hand, only about 25% of the collected PM 2.5 data were below the PM 2.5 annual threshold limit at stations LB, BM and MM.This indicates that high PM 2.5 levels were present in the areas surrounding these three monitoring stations virtually throughout the study period, exposing residents to health risks associated with PM 2.5 .
Similar studies on ambient air quality assessment around brick sites across the world showed moderate to high levels of PM 2.5 at locations close to the brick kilns, which in most cases exceeded the NAAQS permissible limit of PM 2.5 of the respective countries.Saha et al. 92 in their study involving ambient PM 2.5 monitoring around 12 selected brick kiln clusters in Rajashahi and Gazipur districts, between January to April, 2016 and 2017, obtained mean PM 2.5 concentration range of 25.7-2298 µgm −3 across the selected sites.Similarly, high PM 2.5 concentration range between 190 and 2305 µgm −3 was reported by Subhanullah et al. 93 in a 2-days ambient monitoring of PM 2.5 around brick kilns in Northern Pakistan.PM 2.5 assessment studies around the brick kilns in Asia were not conducted within a year cycle, but only considered the active periods and operational hours of the brickmaking industries, hence, higher level of PM 2.5 due to brickmaking were reported 49,92,93 .In addition, majority of the brick kiln used in those studies are Fixed Chimney Brick Kilns which are associated with relatively high release of PM 2.5 compared to the South African clamp kiln 49,50 .

Diurnal variability of PM 2.5
PM 2.5 concentrations from the monitoring stations situated in the two groups showed similar diurnal trends (Fig. 5a).Two PM 2.5 peaks were observed; in the morning (6h00-9h00), and evening (18h00-21h00) throughout the study period.The observed morning and evening PM 2.5 peaks coincide with periods of high traffic congestion, biomass combustion for domestic and industrial purposes such as brick firing.Furthermore, the stable atmospheric condition in the evening (18h00) over the early hours of the day (6h00) could also have aggravated the extent of pollution around the brick kilns in Vhembe District.Similarity in the PM 2.5 patterns generated by the upwind and downwind groups indicates the possible release of PM 2.5 from similar sources.According to Vhembe District Municipality 94 (2020/21) report, fugitive emissions from the brickmaking industries especially the clamp kiln is a major source of PM 2.5 in Vhembe District, during winter.Other primary sources of PM 2.5 emissions include the commercial and small-scale industries, emissions from vehicles, domestic fuel combustion, dust entrainment, biomass combustion, waste disposal and wildfire 94 .The higher peak observed in the evening shows possible increased frequency of activities from PM 2.5 sources at night.From the survey and interviews conducted with the brick kiln workers, brick firing was claimed to be carried out by majority of the brickmakers mostly in the evening when the ambient temperature is low, and occasionally in the morning 51 .Regarding the PM 2.5 levels in monitoring stations in the upwind and downwind direction, concentrations of PM 2.5 at the downwind direction for some of the stations were lower than monitoring station (DR and MG) situated in the upwind direction.This gives an indication of presence of more prominent www.nature.com/scientificreports/PM 2.5 source(s) upwind the monitoring stations DR and MG which could be responsible for the relatively higher PM 2.5 concentrations.

Diurnal variability of PM 2.5 across the seasons
To gain a better understanding on the PM 2.5 trends, their temporal variation across the communities of the brick kilns, the compliance levels as well as the potential sources, it would be needful to sub-categorise the diurnal plot based on seasons.The seasonal diurnal variabilities of PM 2.5 with respect to their compliance levels to both the daily and annual NAAQS threshold of 20 and 40 µgm −3 , respectively, are presented in Fig. 5b.Generally, two PM 2.5 concentration peaks (morning and evening peaks) were observed across all seasons at all the monitoring stations with the evening peaks being more prominent.It is noteworthy that the morning peaks of PM 2.5 at all the monitoring stations fell below the daily NAAQS all through the study period, except for few (three) monitoring stations during winter.
Based on the compliance levels of PM 2.5 at the monitoring stations during autumn (March-May 2021), PM 2.5 concentrations at all the monitoring stations fell below the annual and daily NAAQS thresholds between 13h00 and 15h00.Additionally, all the monitoring stations were in compliance with the daily NAAQS for PM 2.5 between 0h00 and 17h00, while the evening peak in majority of the monitoring stations was non-compliant with both the annual and daily NAAQS threshold limits.In winter (June-August 2021), the highest level of diurnal PM 2.5 concentrations were recorded at all the monitoring stations.In addition, the PM 2.5 concentrations at all monitoring stations exceeded the annual NAAQS during the day except for MG, which fell below it between 12h00 and 15h30.The evening PM 2.5 peaks (17h00-0h00) observed at all the monitoring stations exceeded the daily NAAQS limit.
During spring (September-November 2021), the intensities of some of the observed prominent sources of PM 2.5 in the study area dropped drastically.For instance, the industrial operation of the brick kilns was greatly reduced in spring due to unfavourable weather conditions.Also, the combustion of biomass for space heating purpose is brought to a halt.Unlike the other seasons where the morning and evening peaks for all the monitoring stations occurred about the same period, varying PM 2.5 peaks were observed at all the stations during the morning and the evening, with the variation in the evening peaks being more noticeable in spring.The drastic reduction in some of the prominent PM 2.5 sources such as brickmaking, biomass combustion for space heating, might be responsible for varying PM 2.5 peaks observed across the monitoring stations.Similar evening peaks (15h00-18h00) and (19h00-21h00) were observed at monitoring stations (MR and TP) and (MM, DR and LB), respectively.MR and TP are educational parastatals, high levels of PM 2.5 in these monitoring stations correlates with peak periods of movement of vehicle conveying the student and staff members from school.On the other hand, MM and LB are situated close to T-junctions linked to network of roads.Maximal level of PM 2.5 between 19h00 and 21h00 signifies peak period of vehicular congestion.DR is a lodge, peak PM 2.5 concentration between 19h00 and 21h00 could be possibly from PM 2.5 emission from vehicles bringing in the clients as this period corresponds with the arrival time of most of the customers.During summer, generally low levels of below NAAQS daily threshold limit of 40 µgm −3 were reported in all the monitoring stations.However, from 6h00 to 15h00, majority of the monitoring stations fell below the annual NAAQS threshold limit.The NAAQS annual limit for PM 2.5 was exceeded at all the monitoring stations only during the morning and evening peaks of PM  www.nature.com/scientificreports/one/ two PM 2.5 peaks, coinciding with peak periods of PM 2.5 generating activities [95][96][97][98][99][100][101] .The summary of PM 2.5 diurnal pattern in similar studies are presented in Table 6.

Meteorological conditions and bivariate polar plot
Meteorology has a great influence on PM 2.5 concentration, and its transportation from the source(s) to other locations in Vhembe District.Wind speed and direction provides directional information on the source(s) of PM 2.5 .Figure 6a(i) and (ii) illustrate the annual and seasonal wind roses for the study area during the period of study.Throughout the period of study, inflow of wind was predominantly from the south-east with the maximum wind speed from this direction in summer.High wind speed brings about the dispersion of PM 2.5 away from the source(s), thus leading to the dilution and reduction in PM 2.5 concentration within the given locality.On the other hand, reduced wind speed in a specified direction results in the accumulation of PM 2.5 close to the source.The influence of wind speed and direction on PM 2.5 concentrations across the monitoring stations were assessed using the bivariate polar plots.The polar plot shows quite similar patterns for all the monitoring stations.Overall, high PM 2.5 concentrations were observed at all the stations at the center, North, west and southwest, although at varying concentrations across the monitoring stations (Fig. 6b(i) & (ii)).At low to moderate wind speed of 1-5 ms −1 , generally low PM 2.5 levels ≤ 25 µgm −3 were observed at monitoring stations MR, MT, TP, moderate levels < 30, at station ST, and relatively high PM 2.5 levels > 30 µgm −3 at stations MG, DR, MM, LB and BM in the aforementioned directions compared to the other regions around the monitoring stations.The high PM 2.5 at the center, the northern, western and southwestern areas of the monitoring stations indicate the potential influence of the weak wind speed on PM 2.5 concentration in these directions.Low wind speed facilitates the accumulation of PM 2.5 released from local sources around the monitoring stations.The high concentration experienced close to the center at low wind speed ≤ 2 ms −1 at all the monitoring stations indicate PM 2.5 levels from dominant local sources such as vehicular emission, and domestic biomass combustion.For all the stations, R534 arterial road situated northerly is a potential source while the high concentrations observed at south-western direction at wind speed > 6 ms −1 indicate possible release of PM 2.5 from high rising stack of industries located in this direction.
For monitoring stations (MG and DR) situated upwind the brick kilns, similar trends in the polar plots were observed; highest PM 2.5 concentrations > 35 µgm −3 were observed at the center and in the west direction at low to moderate wind speed of 1-5 ms −1 .Potential PM 2.5 sources at the center include domestic biomass and waste burning, bush burning, and vehicular source.Conversely, Thohoyandou township and Shayandima industrial estate situated in the west and southwest, respectively, are potential PM 2.5 sources responsible for high PM 2.5 levels at higher wind speed (> 6 ms −1 ) at these monitoring stations.According to the source apportionment study on PM 2.5 conducted by Novela et al. 82 , urban emission and industrial sources account for 9.56% and 21.06%, respectively, of PM 2.5 generated in Vhembe District.
For monitoring stations, downwind, the brick kiln, PM 2.5 levels at monitoring stations MM, LB and BM showed similar trends.High concentrations > 30 µgm −3 were observed at the center, the northern and western regions of these monitoring stations.The high PM 2.5 concentrations displayed at the center of monitoring stations MM and LB could possibly be from networks of road surrounding the monitoring stations.Recent study has reported about 10.31% of PM 2.5 emissions originating from vehicular emissions within Vhembe District 82 .Conversely, the high PM 2.5 concentration observed at station BM could possibly be attributed to domestic fuel combustion.Additionally, there were clusters of informal brick making industries which employ the use of woodbased clamp kiln for brick firing situated westly and south-westly the monitoring station BM.The clusters of informal brickmaking industries are potential sources responsible for the high PM 2.5 in west and southwestern regions of station BM.Previous studies have revealed that wood combustion particles accounts for about 24.5% of particulate matters in Limpopo 102 .For monitoring stations LB and ST, construction sites situated at the southwestern region of these monitoring stations are potential sources, however, PM 2.5 levels are mainly localised at monitoring stations MR, MT and TP.

Limitation of the study
In this study, the regulatory standard required for quality assurance of PM 2.5 data was lowered to accommodate the modest functionality of Dylos DC1700.For the collating PM 2.5 data, 70% completion was used for the daily and monthly data collation instead of the approved data completion level of 80%.Whereas, for collation of the seasonal data, collated season must have at least one valid month.Invariably, all the monitoring stations were not represented throughout the seasons in the study period due to limited available data.
Other limitations include: • Comparison of the monitoring sites was done based on unequal levels of data availability.
• Collocation of the Dylos device was only conducted at the North-west university.
• The collocation exercise was conducted for 5 months (June-October 2021) covering only two seasons (winter and spring) of the year.• The PM 2.5 data collected from the Dylos devices were not corrected to account for relative humidity (RH).In the future, it would be beneficial to investigate the impact of RH on the correction factor of Dylos DC1700.

Conclusion
PM 2.5 is of great challenge to scientists all over the world due to its short-and long-term effects on human health.Monitoring of PM 2.5 creates awareness to the public about air quality around them, helping them to take necessary precautions.The conventional method involving the use of federal reference or federal equivalent methods is the most acceptable method.However, not all regions of the world can afford it.To enhance the understanding of the air quality around the less privileged regions of the world, the use of network of well calibrated low-cost sensors is a necessity.For the first time in Africa, a correction factor for PM 2.5 readings from Dylos DC1700 particulate matter sensor was generated.Additionally, the long-term applicability of Dylos for continuous monitoring of outdoor PM 2.5 was assessed.Findings from the study have shown that all the PM 2.5 monitoring stations were non-compliant with the annual NAAQS threshold limit of 20 µgm −3 for PM 2.5 .However, about 50% of the PM 2.5 concentrations at majority of the stations fell below the NAAQS annual threshold limit.Although, none of the monitoring stations exceeded the NAAQS thresholds in summer, highest PM 2.5 concentrations and daily exceedances of the daily NAAQS limit were observed during winter across all the monitoring stations.The high PM 2.5 level experienced especially in winter could be drastically reduced through interventions such as electrification, the provision of clean energy alternatives at subsidised rates, and public enlightenment.The polar plots of all the monitoring stations show similar trends with high PM 2.5 concentrations majorly at the centre, north, west and south-west areas of the monitoring stations.PM 2.5 concentrations were majorly localised at monitoring stations MR, MT and TP.The polar plot also showed maximal impact on PM 2.5 concentration from brick kilns on monitoring station BM located downwind and ≤ 500 m from the brick kilns.For the other monitoring stations, influx of PM 2.5 were observed at wind speed > 5 ms −1 at the southwest, signifying industries as the potential PM 2.5 source.Other likely potential sources responsible for the high PM 2.5 concentrations around the monitoring stations include construction site, major R534 road and Thohoyandou township.The reduced impact of PM 2.5 emission from brickmaking industries at most of the monitoring stations displayed by the polar plot does not necessarily means PM 2.5 emission from the brickmaking industries has no influence on the ambient PM 2.5 level, however, there is high possibility of increased PM 2.5 readings in areas very close to the brick kilns which were not included as monitoring stations.In conclusion, the application of Dylos DC1700 offers an excellent solution to bridging the data gap and establishing air pollution control policies for public health protection in areas with sparse air quality data.

Figure 1 .
Figure 1.The study area map showing the geographical locations of the monitoring stations produced using ArcGIS 10.8.

Figure 2 .
Figure 2. The monitoring equipment (a) Modified air monitoring device.(b) Installed air monitoring device connected to a laptop.

Figure 3 .
Figure 3. Relationship between the GRIMM EDM 180 and Dylos DC1700 readings during the collocation exercise.

Figure 4 .
Figure 4. Temporal variation in PM 2.5 .(a) Seasonal variation in PM 2.5 concentration of the monitoring stations in the upwind direction of the brickmaking industries and the standard deviations (bars).(b) Seasonal variation in PM 2.5 concentration of the monitoring stations in the downwind direction of the brickmaking industries and the standard deviations (bars).(c) Annual ambient PM 2.5 compliance status of monitoring stations close to the brick kilns in Vhembe District.The box-and-whisker plots show the minimum and maximum values (the lower and higher end of the whiskers); means (red dots); 25th, 50th, and 75th quartiles; and the annual NAAQS threshold of PM 2.5 (solid red line).

Figure 5 .
Figure 5.The diurnal variation of PM 2.5 across the monitoring stations (a) Annual variability (b) Seasonal variability.

Received: 15 Figure 6 .
Figure 6.Meteorological conditions and bivariate polar plot.(a)(i) Annual wind rose of hourly data in Thulamela municipality.(a)(ii) Seasonal windrose of hourly data in Thulamela municipality (March 2021 to February 2022) indicating the prevailing wind direction, the mean windspeed and the percentage calmness.(b) (i) Polar plots of mean concentrations of PM 2.5 for the monitoring period at the monitoring stations upwind the brickmaking industries.(b)(ii) Polar plots of mean concentrations of PM 2.5 for the monitoring period at monitoring stations downwind the brickmaking industries.

Table 1 .
The selected sites and their characteristics.

Table 2 .
Data completeness for all the PM 2.5 monitoring station from March 2021-February 2022.Data availability with percentages less than 70% are emboldened.

Table 3 .
82scriptive statistics for the daily averaged ambient PM 2.5 mass concentrations in µgm-3 at nine monitoring stations in Vhembe District, for the study period March 2021 to February 2022, sub-categorised by season.The rows of data italicised are the seasons.-industrialisedsetting,that of Novela et al.82was in a control environment with limited PM 2.5 emission sources.This study as well as other similar studies enlisted in Table5indicated reduced levels of PM 2.5 during rainy season (spring and summer) and relatively higher PM 2.5 concentrations during dry season (autumn and winter).

Table 4 .
Variation in PM 2.5 across the sampling months (mean ± sd) in microgram per cubic meter.

Table 5 .
Comparison of result from Vhembe District to that from other cities in Africa.

Table 6 .
Comparison of PM 2.5 diurnal variability pattern in Vhembe District with similar studies.